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Abstract 

This paper studies the problem of option replication in general stochastic volatility markets 
with transaction costs, using a new specification for the volatility adjustment in Leland’s 
algorithm [23] . We prove several limit theorems for the normalized replication error of Leland’s 
strategy, as well as that of the strategy suggested by Lepinette [27]. The asymptotic results 
obtained not only generalize the existing results, but also enable us to fix the under-hedging 
property pointed out by Kabanov and Safarian in [18]. We also discuss possible methods to 
improve the convergence rate and to reduce the option price inclusive of transaction costs. 
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1 Introduction 

Leland [23] suggests a simple method for pricing standard European options in markets with 
proportional transaction costs. He argues that transaction costs can be accounted for in the 
option price by increasing the volatility parameter in the classical Black-Scholes model [4]. Leland 
then claims, without giving a mathematically rigorous proof, that the replicating portfolio of the 
corresponding discrete delta strategy converges to the option payoff as the number of revisions n 
goes to infinity, if the transaction cost rate is a constant independent of n, or decreases to zero at 
the rate . The latter statement is proved by Lott in his PhD thesis [30]. In fact, this property 
still holds if the transaction cost coefficient converges to zero at any power rate [18]. 

However, a careful analysis shows that the replicating portfolio does not converge to the option 
payoff when the cost rate is a constant independent of n. Kabanov and Safarian [18] find an 
explicit limit for the hedging error, which is negative, showing that the replication problem is 
not completely solved in Leland’s framework. Pergamenshchikov [34] obtains a weak convergence 
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for the normalized hedging error and points out that, for the case of constant transaction cost, 
the rate of convergence in Kabanov-Safarian’s result is This limit theorem is of practical 

importance because it provides the asymptotic distribution of the hedging error. Note that the 
rate of convergence can be improved using non-uniform revisions [27, 9]. In these papers, Lepinette 
and his co-authors suggest a modification to Leland’s strategy to solve the discrepancy identified 
by Kabanov and Safarian. For a recent account of the theory, we refer the reader to Section 2 and 
[18, 24, 26, 27, 13, 14, 9, 34], 

In this study, we examine the problem of approximate hedging of European style options in 
stochastic volatility (SV) markets with constant transaction costs (the reader is referred to e.g. 
[11] and the references therein for motivations and detailed discussions related to SV models). In 
particular, we establish a weak convergence for the normalized hedging error of Leland’s strategy 
using a simple volatility adjustment, in a general SV setting. The results obtained not only 
generalize the existing results, but also provide a method for improving the rate of convergence. 
Furthermore, it turns out that superhedging can be attained by controlling a model parameter. 

Let us emphasize that the classic form for adjusted volatility proposed in [23] and applied in 
[18, 19, 24, 25, 27] may not be applicable in SV models. The reason is that option pricing and 
hedging are intrinsically different in SV markets than in the classical Black-Scholes framework. 
In particular, the option price now depends on future realizations of the volatility process. In 
general, this information may not be statistically available for all investors. To treat this issue, we 
suggest a new specification for adjusted volatility in Leland’s algorithm. Although we employ an 
artificially modified volatility, simpler than the well-known version used in the previous literature, 
the same asymptotic results are obtained for SV contexts. In addition, the rate of convergence 
can be improved by controlling a model parameter. Note that, in the above-mentioned papers, 
approximation procedures are mainly based on moment estimates. This essential technique no 
longer works in general SV models, unless some intrinsic conditions are imposed on the model 
parameters [2, 28]. It is useful to remember that our goal is to establish a weak convergence for the 
normalized replicating error which only requires convergence in probability of the approximation 
terms. Thus, in the approximation procedure, the integrability issue can be avoided in order to 
keep our model setting as general as possible. 

As discussed in [34], the option price (inclusive of transaction costs) in Leland’s algorithm 
may be high (it, in fact, approaches the buy-and-hold price), even for small values of the revision 
number. Another practical advantage of our method is that the option price can be reduced as 
long as the option seller is willing to take a risk in option replication. This approach is inspired by 
the theory of quantile hedging [10]. 

The remainder of the paper is organized as follows. In Section 2, we give a brief review of 
Leland’s approach. Section 3 is devoted to formulating the problem and presenting our main 
results. Section 4 presents some direct applications to pricing and hedging. Section 5 discusses 
common SV models that fulfill our condition on the volatility function. A numerical result for Hull- 
White’s model is also provided for illustration. Section 6 connects our results to high-frequency 
markets with proportional transaction costs. The proofs of our main results are reported in Section 
7. Auxiliary lemmas can be found in the Appendix. 


2 Approximate hedging with transaction costs: A review of 
Leland’s approach 

In a complete no-arbitrage model (i.e., there exists a unique equivalent martingale measure under 
which the stock price is a martingale), options can be completely replicated by a self-financing 
trading strategy. The option price, defined as the replication cost, is the initial capital that the 
investor must invest to obtain a complete hedge. In fact, the option price can be computed as 
the expectation of the discounted claim under the unique equivalent martingale measure. This 
principle plays a central role in the well-known Black-Scholes model. For simplicity, let us consider 
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a continuous time model of a two-asset financial market on the time interval [ 0 , 1 ], where the bond 
price is equal to 1 at all times. The stock price dynamics follow the stochastic differential equation 

dSt = a^S^dW ^, S'o given, ( 1 ) 

where S'o and (Jg are positive constants and is a standard Wiener process. As usual, 

let J't = a{W^ , 0 < u < t}. We recall that a financial strategy (/3j,is an admissible 
self-financing strategy if it is bounded from below, (Jy) - adapted with j^{\Pt\ + 7 t)dt < oo a.s., 
and the portfolio value satishes 

= /?* + ItSt = ^0+ [ 7 „dS„, t G [0,1]. 

do 

The classic hedging problem is to find an admissible self-financing strategy (/ 3 (, 7 () whose terminal 
portfolio value exceeds the payoff h{Si) = max(S;^ — K, 0), or 

Vi = Vo+ [ 7 ndS„ > h(Si) a.s., 
do 

where K is the strike price. The standard pricing principle shows that the option price C(t,S^) is 
given by the well-known formula [4] 

C{t,x) = C(t,x,aQ) = x^{v{t,x)) — K^(y{t,x) — crg-\/l — t ), (2) 

where 

v{t,x) = v{al{l-t),x) and v(A,a;) = ^ ■ (3) 

Here, $ is the standard normal distribution function. In the following, we denote by (p the Af{0, 1) 
density: (p{z) = One can check directly that 

C^{t,x) = <i>(y{t,x)) and C^^{t,x) = ( 4 ) 

XCTgVl — t 

By assuming that continuous portfolio adjustments are possible with zero transaction costs. Black 
and Scholes [4] argue that the option payoff can be dynamically replicated using the delta strategy 
(i.e., the partial derivative of the option price with respect to the stock price). 

It is clear that the assumption of continuous portfolio revision is not realistic. Moreover, 
continuous trading would be ruinously expensive in the case of nonzero constant proportional 
transaction costs because the delta strategy has infinite variation. This simple intuition contradicts 
the argument of Black and Scholes that, if trading takes places reasonably frequently, then hedging 
errors are relatively small. Therefore, option pricing and replication with nonzero trading costs 
are intrinsically different from those in the Black-Scholes setting. Note that it may be very costly 
to assure a given degree of accuracy in replication with transaction costs. In what follow, we show 
that Leland’s increasing volatility principle [23] is practically helpful in such contexts. 

2.1 Constant volatility case 

Leland’s approach [23] provides an efficient technique to deal with transaction costs. This method 
is simply based on the intuition that transaction costs can be accounted for in the option price 
as a reasonable extra fee, necessary for the option seller to cover the option return. It means 
that in the presence of transaction costs, the option becomes more expensive than in the classic 
Black-Scholes framework. This is intuitively equivalent to an increase in the volatility parameter 
in the Black-Scholes formula. Let us shortly describe Leland’s approach [23, 18]. Suppose that for 
each trading activity, the investor has to pay a fee directly proportional to the trading volume. 
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measured in dollar value. Assume that the transaction cost rate is given by the law where 

n is the number of revisions. Here, 0 < a < 1/2 and > 0 are two fixed parameters. The basic 
idea of Leland’s method is to replace the true volatility parameter in the Black-Scholes model by 
a, artificially modified as 

= cTp + with Q = . (5) 


In this case, the option price is given by C{t,x) = C{t,x,a), the Black-Scholes’s formula. For the 
problem of option replication, Leland suggests the following discrete strategy, known as Leland’s 
strategy. 




2=1 


ti = ie {l,2,..,n}. 
n 


( 6 ) 


Here, the number of shares held in the interval is the delta strategy calculated at the left 

bound of this interval. Then, the replicating portfolio value takes the form 


y" = + 



— K^n 




(7) 


where the total trading volume is ~ It I (measured in dollar value). Recall 

that the option price C{t,x) is the solution of the Black-Scholes PDF with the adjusted volatility 
a 

+ ^a‘^x'^d^^{t,x) =0, 0<t<l; d{l,x) = h{x). (8) 

Using Ito’s formula, we can represent the tracking error, U" — h{Si), as 



C,{t,SMdS, 






(9) 


Remark 1 (Leland). The specific form (5) results from the following intuition: the Lehesgue inte¬ 
gral in (9) is clearly well approximated by the Riemann sum of the terms ^)At, 

while 4 1 7 " ~ 7 " I replaced by 




^xa;(^i-l)~ 'XQ\/2l{n'K) ^ 


because E|AH/ | = = \f2jijrri). Hence, it is reasonable to expect that the modified 

volatility defined in (5) will give an appropriate approximation to compensate transaction costs. 

Leland [23] conjectures that the replication error converges in probability to zero as n —> 00 for 
the case of constant proportional transaction cost (i.e., a = 0). He also suggests, without giving 
a rigorous proof, that this property is also true for the case a = 1/2. In fact, Leland’s second 
conjecture for a = 1/2 is correct and is proved by Lott in his PhD thesis [30]. 

Theorem 2.1 (Leland-Lott [23, 30]). Fora = 1/2, strategy (6) defines an approximately replicat¬ 
ing portfolio as the number of revision intervals n tends to infinity 


P- lim y|" = /i(4). 

n—foo 

This result is then extended by Ahn et al. in [1] to general diffusion models. Kabanov and Safarian 
[18] observe that the Leland-Lott theorem remains true as long as the cost rate converges to zero 
as n —>■ 00 . 

Theorem 2.2 (Kabanov-Safarian [18]). For any 0 < a < 1/2, P — lim„_,,^ U" = h{Si). 
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In [26, 19], the authors study the Leland-Lott approximation in the sense of convergence 
for the case a = l/ 2 d 

Theorem 2.3 (Kabanov-Lepinette [26]). Let a = 1/2. The mean-square approximation error of 
Leland’s strategy, with g defined in (5), satisfies the following asymptotic equality 

El {V^ — h{Si)Y = An~^o{n~^) n —>■ oo, 

where A is some positive function. 

Theorem 2.3 suggests that the normalized replication error converges in law as n —)■ oo. 

Theorem 2.4 (Lepinette-Kabanov [19]). For a = 1/2, the processes T" = — h{Si)) 

converge weakly in the Skorokhod space I2[0,1] to the distribution of the process Y, = fj B(Sf.)dZ^, 
where Z is an independent Wiener process. 

Remark 2. An interesting connection between this case and the problem of hedging under propor¬ 
tional transaction costs in high-frequency markets is discussed in Section 6. 

It is important to note that Leland’s approximation in Remark 1 is not mathematically correct 
and thus, his first conjecture is not valid for the case of constant transaction costs. In fact, as 
n ^ oo, the trading volume can be approximated by the following sum (which converges in 
probability to J(S'j, p) defined in ( 11 )) 

n 

- E TiK-1, S,^J\a,Q-^Z, + q{\_„S,^J\ AA,, 

where = A^. = d‘^{l - tf), Z^ = and 

x) = V3 (v(A, a;)), q{X, x) = ^ ■ (10) 

A careful study confirms that there is a non trivial discrepancy between the limit of the replicating 
portfolio and the payoff for the case a = 0 . 

Theorem 2.5 (Kabanov-Safarian [18]). For a = 0, R" converges to h{Si) + min(S']^,Ar) — 

K„,J{Si,g) in probability, where 

p + (X> 

J{x,g)=x / X~^/‘^ip{\,x)E \gZ-\-q{\,x)\ dX, (11) 

•to 

with g = crQg~^ and Z ~ Af{0, 1) independent of Si. 

Under-hedging: It is important to observe that the problem of option replication is not com¬ 
pletely solved in the case of constant transaction costs. Indeed, considering that E \gZ\ = 1/{2 k.^) 
and the identity 

OO 

X j X~^/^ ip{X,x)dX = 2mh\{x,K), (12) 

0 

we obtain (for the parameter g given in (5)) that min(a:, K)—k^J{x, g) = equals X~^/'^tp{X, x) (E IgZj — E jgZ -t- 
Now, Anderson’s inequality (see, for example [17], page 155) implies directly that for any q £ R, 


^ Seemingly, mean-square replication may not contain much useful information because gains and losses have 
different meaning in practice. Clearly, if a = 1/2 the modified volatility is independent of n. 
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E \qZ + ( 7 I > E \qZ\ . Therefore, P — (E" — h{Si)) < 0, thus, the option is asymptotically 

under-hedged in this case. 

In approximation procedures, one should also pay attention to the fact that C and its derivatives 
depend on the number of revisions when 0 < a < 1/2. In addition, the coefficient g appearing in 
(5) can be chosen in an arbitrary way. Pergamenshchikov [34] shows that the rate of convergence in 
Kabanov-Safarian’s theorem is n}/'^ and provides a weak convergence for the normalized replication 
error. 


Theorem 2.6 (Pergamenshchikov [34]). Consider Leland’s strategy ( 6 ) with a = 0, and let g in 
(5) be some fixed positive constant. Then, the sequence of random variables 

„i/ 4 (yn _ + n,J{S„g)) (13) 

weakly converges to a centered mixed Gaussian variable as n —>■ 00 . 

Theorem 2.6 is of practical importance because it not only gives the asymptotic information 
about the hedging error, but also provides a reasonable way to fix the under-hedging issue (see 
Section 4) . Darses and Lepinette [9] modify Leland’s strategy in order to improve the convergence 
rate in Theorem 2.6 by applying a non-uniform revision policy (ti)i<i<„, dehned by 

tj = g (i/n) , g{t) = 1 — {1 — t)'^ for some pL > 1. (14) 

The adjusted volatility is then taken as ct/ = (Tp-|-At„crg-\/ 8 / 7 r-^n/'(t), where / is the inverse function 
of g. Furthermore, the discrepancy in Theorems 2.5 and 2.6 can be removed by employing the 
following modified strategy, known as Lepinette’s strategy [27], 


7r 


2 = 1 




(15) 


Theorem 2.7. Let V/" be the terminal value of the strategy (15) with a = 0. Then, for any 
1 < p, < /imax; ^^6 sequence — h{Sfij) weakly converges to a centered mixed Gaussian 

variable as n —)■ 00 , where 




2 (^ + 1 ) 


and /r^ax 


8 


(16) 
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2.2 Time-dependent volatility case 

Assume now that cr is a positive non-random function and the payoff iJ is a continuous function 
with continuous derivatives, except at a finite number of points. Under the non-uniform rebalancing 
plan (14), the enlarged volatility should take the form 

af = a^{t) + K*cr(t)n^/^““\//'(t)8/7r. (17) 

Theorem 2.8 (Lepinette [24]). Let a he a strictly positive Lipschitz and bounded function. More¬ 
over, suppose that H{-) is a piecewise twice differentiable function and there exist a;* > 0 and 
6 > 3/2, such that sup^^^ x^\H (cc)] < oo. Then, for any 1/2 > a > 0, the replicating portfolio 
of Leland’s strategy converges in probability to the payoff H(Si) as n ^ oo. Moreover, for a = 0, 

P - hm Vf = H{S,) + H,{S,) - nMS,), 

n—>-oo 

where H^{-) and i 72 (’) positive functions that depend on the payoff H. 

Remark 3. Theorem 2.7 still holds in the context of Theorem 2.8 (see [27]). 

2.3 Discussion 

From Remark 1, the modified volatility defined by (5) would seem to give an appropriate approxi¬ 
mation that accounts for transaction costs. However, this is not always the case because the option 
price inclusive of transaction costs now depends on the rebalancing number. In more general mod¬ 
els, this specific choice may generate technical issues. For example, in local volatility models [24], 
proving the existence of the solution to ( 8 ) requires patience and effort, because t? depends on the 
stock price. On the other hand, it is interesting to observe that the true volatility plays no role 
in the approximation procedure from a mathematical point of view. In fact, all the results for the 
case 0 = 0 can be obtained by using the form a] = ^/J^t)8>J^:, where the first term 

cr^(t) has been removed. More generally, we can take the following form 

= Q\/nf’{t), (18) 

for some positive constant g, which will be specified later. Of course, the limiting value of trans¬ 
action costs will change slightly. Let us emphasize that using the simple form (18) is important 
for two reasons. First, it allows us to carry out a far simpler approximation than is used in the 
existing literature. Second, Leland’s strategy with a defined in (5) may no longer work in stochas¬ 
tic volatility (SV) markets. Indeed, in those markets, option prices depend on future volatility 
realizations, which are not statistically available. We show in the remainder of the paper, that 
the simple form (18) (a deterministic function of t) is helpful for approximate hedging in a very 
general SV setting. It should be noted that the approximation methodology developed here still 
works well for the classical form (5), if the volatility risk premium depends only on the current 
value of the volatility process [36, 37]. 

We conclude the section by mentioning that Leland’s algorithm is of practical importance due 
to its ease of implementation. The case of constant transaction costs a = 0 should be investigated 
in more general situations, for instance, where volatility depends on external random factors, or 
jumps in stock prices are considered. 


3 Model and main results 

Let (H, (Jq)g<;(<j^, P) be a standard filtered probability space with two standard independent 

(•^t)o<t<i adapted Wiener processes (W^^^^) and taking their values in K. Our financial 

market consists of one risky asset governed by the following equations on the time interval [0,1] 

d^, = a(y,)^,dW,(i); dy, = F,{t,y,)dt + ^^(t, y,)(rdIF/i) + Vl-r^dWf)), (19) 
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where — 1 < r < 1 is the correlation coefficient. It is well known in the literature of SDEs that if 


Fi{t,y) and F 2 {t,y) are measurable in {t,y) S [0,r] x M, linearly bounded and locally Lipschitz, 
there exists a unique solution y to the last equation of system (19). For this fundamental result, 
see Theorem 5.1 and [12, 29]. For simplicity, assume that the interest rate equals zero. Thus , the 
non-risky asset is chosen as the numeraire. 

In this section, we consider the problem of approximate hedging with constant proportional 
costs using the principle of increasing volatility for model (19). As discussed in Subsection 2.3, the 
adjusted volatility is chosen as 



( 20 ) 


The replicating portfolio is revised at (t^), as defined by (14). The parameter p > 0 plays an 
important role in controlling the rate of convergence and is specified later. As shown below, the 
limiting value of the total trading volume is essentially related to the dependence of g on the 
number of revisions. 

Remark 4. Intuitively, using an independent adjusted volatility seems unnatural because it fails to 
account for market information. However, the techniques developed in this note are well adapted 
to the case where the adjusted volatility depends on a volatility process driven by an independent 
Brownian motion. In such a context, if the volatility risk premium depends only on the current 
volatility process, then the no-arbitrage option price (without transaction costs) is the average of 
the Black-Scholes prices over the future paths of the volatility process [36, 37]. 

Recall that C{t,x) is the solution of the Cauchy problem ( 8 ) with two first derivatives, as 
described in (4): C„.{t, x) = $(v(Aj, x)) and C^,,.(t, x) = bp (A^, x), where 



( 21 ) 


t 


Remark 5. Section 4 will show that the under-hedging situation pointed out in [18] can be fixed 
by controlling the parameter g. 

We make use of the following condition on the volatility function. 

(C]^) Assume that a is a C'^ function and there exists cr^in such that 


0 < cr^in < cr{y) for ally GR and sup F,[a'^{y^) + \cr'{yt)\] < oo. 


0<t<l 


Assumption (Ci) is not restrictive and is fulfilled in many popular SV models (see Section 5 and 
[35]). 


3.1 Asymptotic results for Leland’s strategy 

Let us study the replication error for Leland’s strategy 7 " defined in ( 6 ). The replicating portfolio 
F" is defined by (7). Now, by Ito’s formula, 



0 


( 22 ) 


where {a[ - a^iyt)) S[d,,,,{t, St)dt. Setting Vq = 0(0, Sq), 


we can represent the repli¬ 


cation error as 



(23) 


where /2 „ = ( 7 " — C,j,{t, dA^ and J„ is dehned as in (7). 
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Let us first emphasize that complete replication in SV models is far from obvious. In our setting, 
I 2 n converges to zero faster than n^, with /3 defined as in (16). The gamma error ^ approaches 
2 minif) at the same rate. On the other hand, the total trading volume converges in 
probability to the random variable g), defined by 

|•+oo 

A"^/^^(A,x)E |cr(j/)p"^Z + 9(A,a:)| dA, (24) 

"'O 

where Z ~ A/'(0,1) independent of Si and yi. Our goal is to study the convergence of the normalized 
replication error corrected by these explicit limiting values, by applying the theory of limit theorems 
for martingales [15]. To do so, we search for the martingale part in the approximation of the above 
terms by developing a special discretization procedure in Section 7. 

Theorem 3.1. Suppose that condition (Ci) holds and g > 0 is a fixed positive constant. Then, 

n^iV” - h{Si) - inm{Si,K) + K,J{Si,yi,g)) 
weakly converges to a centered mixed Gaussian variable as n ^ oo. 

Remark 6. This theorem is a generalization that includes an improved convergence rate of the 
results in [18, SfJ, where the uniform revision is taken and the volatility is assumed to be a constant. 

Remark 7. Note that h{x) + m.in{x,K) = x, where h{x) = {x — K)'^ is the payoff of a classical 
European call option. Then, from Theorem 3.1, the wealth process Vf' approaches Si—K^J{Si,yi, g) 
as n ^ oo. This can be explained by the fact that the option is now sold at a higher price because 
C{0,SQ,a) Sq as a ^ oo. In other words, Leland’s strategy now converges to the well-known 
buy-and-hold strategy [22]: buy a stock share at time t = 0 for price Sq and keep it until expiry. 

We now present a method for improving the rate of convergence in Theorem 3.1. To this end, 
by letting p —)■ oo, we observe that 

|■+ao 

lim J{x,y,g) = x / A“^/^^(A, a:)|g(A, a::)|dA :=J*{x), (25) 

e-s-oo Jq 

which is independent of y. This suggests that the rate of convergence in Theorem 3.1 can be 
improved if g is taken as a function of n. Our next result is established under the following 
condition on g. 

(C 2 ) The parameter g = g(n) is a function of n such that 

_ /i. 

lim g{n) = 00 and lim gn ^((^+ 2 ) =0. 

n—>-oo n—>-oo 

Theorem 3.2. Under conditions (C^^) — (C 2 ), 

0^{V[[^-hiS,)-uim{S„K) + K,J*iS,)), with 9^=n^g^P, 
weakly converges to a centered mixed Gaussian variable as n —>■ 00 . 

Remark 8. The asymptotic distributions in Theorems 3.1 and 3.2 are explicitly determined in the 
proofs in Section 7. 

3.2 Asymptotic result for Lepinette’s strategy 

Let us study the replication error of Lepinette’s strategy 7 ", as defined in (15). As before, the 
replicating portfolio is E” = Eg + f^JfdSf — k^J^, where 

n 

= (26) 
2=1 
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Now, by Ito’s formula, the tracking error is 

( 27 ) 

where / 2 _„ = /q’”^ C'„,j(u, S'^)du. Then, we have the following strengthening 

of Theorem 2.7. 

Theorem 3.3. Suppose that (Ci) is fulfilled. Then, for any g > 0, the sequence 
{V^ — h{Si) — r]min{Si, K)), with rj = I — K.^a{yi)g~^ y/S/n, 
weakly converges to a centered mixed Gaussian variable as n —>■ oo. 

Remark 9. Theorem 2.7 can he established from Theorem 3.3 with g = Ks,<7^f%pK when the 
volatility is a constant. In addition, in our model, the parameter p takes its values in the interval 
[1,2), which is slightly more general than the condition imposed in Theorem 2.1. Moreover, if the 
classical form of adjusted volatility is applied for Lepinette’s strategy, then complete replication can 
he reached by taking g = K^y^S/n, and we again have the result established in [9]. 

Corollary 3.1. Under conditions (Ci) — (C 2 ), the wealth sequence 17^ converges in probability to 
/i(S'i)+min(5i,X) = S'i. 

Note that we do not obtain an improved convergence version of Theorem 3.3 because k* con¬ 
verges to zero at the order of g. 

4 Application to the pricing problem 

In this section, we present an application to the problem of option pricing with transaction costs. 
We first emphasize that it is impossible to obtain a non-trivial perfect hedge in the presence of 
transaction costs, even in constant volatility models. In fact, the seller can take the buy-and-hold 
strategy, but this leads to a high option price. We show below that the price can be reduced in 
certain ways so that the payoff is covered with a given probability. 

4.1 Super-hedging with transaction costs 

To be on the safe side, the investor searches for strategies with terminal values greater than the 
payoff. Such strategies are solutions to dynamic optimization problems. More precisely, let H he a. 
general contingent claim and let A{x) and be the set of all admissible strategies tt with initial 
capital X and the terminal value of strategy tt, respectively. Then, the super-replication cost for 
H is determined as 

C/q = inf {x G K : d-TT e A{x), > H a.s.} , (28) 

(see [22] and the references therein for more details). In the presence of transaction costs, Cvitanic 
and Karatzas [8] show that the buy-and-hold strategy is the unique choice for super-replication, 
and then Sq is the super-replication price. In this section, we show that this property still holds 
for approximate super-hedging. The following observation is a direct consequence of Theorem 3.2 
when p is a function of n. 

Proposition 4.1. Under conditions (Ci) — (C 2 ), P — lim„_,,,,^ > h{Si). The same property 

holds for Lepinette’s strategy. 

Proof. Note first that J*{x) < mm{x, K), for all x > 0. Hence, by Theorem 3.2, 

P - lim {Vf- - h{Si)) > (1 - K,)mm{Si,K). (29) 

n—>-oo 

The left-hand side is obviously non-negative as < 1. The conclusion follows from Theorem 
3.3. □ 
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4.2 Asymptotic quantile pricing 

As seen ealier, super-hedging in the presence of transaction costs leads to a high option price. 
Practically, one can ask by how much the initial capital can be reduced in exchange for a shortfall 
probability at the terminal moment. More precisely, for a given significance level 0 < e < 1, the 
seller may look for hedges with a minimal initial cost 

inf {x G M, dTT e A{x) : P >H)>l-e}. 

This construction is motivated by quantile hedging theory, which goes back to [10, 33]. For related 
discussions, we refer to [10, 33, 34, 5, 7, 6]. Here, we adapt this idea to the hedging problem. Recall 
that the super-hedging price of Leland’s algorithm is Sq. On the seller’s side, we propose a price 
5Sq < Sq for the option, for a properly chosen 0 < <5 < 1. We then follow Leland’s strategy for 
replication. To be safe at the terminal moment, we need to choose g such that the probability of 
the terminal portfolio exceeding the sum of the real objective (i.e., the payoff) and the additional 
amount (1 — 5)3^ is greater than 1 — e. Here, £ is a significance level predetermined by the seller. 
By Proposition 4.1, this goal can be achieved for sufficiently large g. To determine the option 
price, it now remains to choose the smallest value of S. Motivated by (29), we define this by 

= inf {o > 0 : T(a) > 1 — e} , T(a) = P ((1 — k^) min)^]^, AT) > (1 — 0 ) 50 ). (30) 

Thus, the reduction in the option price is given by (1 — S^)So. Clearly, smaller values of 6^ yield 
cheaper options. 

Next, we show that the option price is significantly reduced, compared with powers of the 
parameter e. 

Proposition 4.2. Assume that = sup^^^ cr(?/) < 00 . Then, for any r > 0 and 5,, defined by 
(30), 


lim(l —5g)£ '" = - 1 - 00 . (31) 

£->■0 

Proof. We first observe that 0 < < 1 and <5^ tends to 1 as £ —>■ 0. Set b = 1 — Then, for 

sufficiently small £ such that > a > 1 — bK/SQ, one has 

1 - £ > P(&min(S'i, X) > (1 - a)So) = 1 - P{SJSo < (1 - a)/b). 


Therefore, 


£ < P {S,/S, < (1 - a)/b) < P (A, < -zj , (32) 

where X^ = a(yf)dWf^^ and = ln(6/(l — a)) — tT^^^/2. To estimate this probability, we 
note that for any integer m > 1, E(A]^)^™ < tT^^(2TO — 1)!! (see [29, Lemma 4.11, p.l30]). Set 
R(v) = 2va'^ . For any 0 < u < l/2o-^ , 

^ ' max ' max’ 


E = y y E (Ai)"™ < y y < _ 

^ m! ^ ^ m\ ’ ~ 1 


^2m , V 
m=0 m—0 

Therefore, for sufficiently small £ > 0, we have 


-R{v) 


e < P(Ai < -zj = P{-X, > zj < Ee’^^i < — 


R{vy 


Then, 1 — a > 6 e where (.^(w) = -y|ln£(l — i?(£’))| /v + a^^^/2. Letting a —?> <5^, we get 

1 — (5^ > be~''‘:^'"\ which implies (31). □ 

The boundedness of the volatility function is essential for the above comparison result. If we 
wish to relax this assumption, the price reduction will be smaller than that in Proposition 4.2. 
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Proposition 4.3. Suppose that Eexp{acr^(j/s)ds} < oo, for some constant a > 1/2. 
for ra = + l)/ 2 a, 


liminf e (1 — <5^) > 0 . 

E —^0 


Then, 

(33) 


Proof. For any positive constant L we set 


T = = inf 


t>0: 


a‘^{ys)ds > L > A 1, 


which is understood to be the first time that the log-price’s variance passes level L. 
(32), 


(34) 


Then, from 


e < P ^(ct) > Ua, J a'^{ys)ds <L^+'p(^J a'^{ys)ds > , (35) 

where S^ia) = e^o '“s/q — k„)/( 1 — a), and i5e > a > 1 — bK/pQ. Note 

that for any p > 0, the stopped process Xt — ^ martingale and = 1. Therefore, 

the first probability on the right side of (35) can be estimated as 

(«J-^’E£;^’(a) = («J-"Exi e^^o < (^xJ-^’e^^ 

where p = {p^ + p)/2. By the hypothesis and Chebysev’s inequality, we have 

P cr'^iys)ds > with C'„=Eexp|ay' tT^(?/,,)ds|. 

Hence, e < {u,,)~p e^^ + By choosing L = a~^ ln(2C'/e) and letting a —>■ Sg, one deduces 

that for any p > 0 and for some positive constant C^, 


I — Se > C^s"^ where 7 *(p) = (p -b l)/( 2 a) + p ^. 

Note that ra = minp^Q 7 *(p) = 7 *(-\/^). Then, including in the last inequality p = \/2a we 
obtain property (33). □ 

Remark 10. It is clear that < 1 for a > 3/2-1- \/2. The condition used in Proposition 4-3 
holds for such a, when a is linearly bounded and y^ follows an Or stein-Uhlenbeck process (see the 
Appendix C). The same quantile pricing result can be established for the Lepinette strategy. 


5 Examples and numerical results 

In this section, we list some well-known SV models for which condition (Ci) is fulfilled. To this 
end, we need some moment estimates for solutions to general SDEs, 

dyt = Fi{t,yf)dt +F2(t,yf)dZ^, p(0) = po, (36) 

where Z is a standard Wiener process and Fi,F 2 are two smooth functions. We first recall the 
well-known result in SDE theory (see for example [12, Theorem 2.3, p.l07]). 

Theorem 5.1. Suppose that Fi(t,y) and F 2 (t,y) are measurable in it,y) G [0,T] x M, linearly 
bounded and locally Lipschitz. //E |pqP™ < oo for some integer m > 1, then there exists a unique 
solution (pj) to (36) and 

E Ip.p™ < (1 -f E |poP'")e“‘, E sup |p,|2™ < M(1 + E \y,\^^), 

0<s<t 

where a,M are positive constants dependent on t,m. 
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In the context of Theorem 5.1, condition (Ci) holds if the volatility function and its derivative 
have polynomial growth \a{y)\ < C(1 + |y|™), for some positive constant C and to > 1. 

Hull-White models: Assume that y^ follows a geometric Brownian motion 

dS't = iVt + and dy^ = y((adt + bdZ^), (37) 

where (T„iin > 0, a and b are some constants, and .Z is a standard Brownian motion correlated with 
W^.. Put y* = supp^j^j^ Ij/J. Then, by Theorem 5.1, we have 

E(y*)2™<C(l + E|%|2-)<oo, TO>1, 

as long as E|yQp™ < oo. Therefore, condition (Ci) is fulfilled in (37). 

Uniform elliptic volatility models: Suppose that volatility is driven by a mean-reverting 
Orstein-Uhlenbeck process 

dS't = (j/f -b a^^^)S^AW^ and dy^ = {a- hy^)dt + dZ. (38) 

In this case, a{y) = y"^ + cr^in- Thus, condition (Ci) is verified throughout Theorem 5.1. 
Stein-Stein models: Assume that 

= yjy1 + S^dw^ and dy^ = (a - feyjdt -b dZ^. (39) 

We have cr(y) = ^/y^~+~o\~ and condition (Ci) is also verified by Theorem 5.1. 

Heston models: Heston [16] proposes a SV model where volatility is driven by a CIR process, 
which also known as a square root process. This model can be used in our context. Indeed, assume 
now that the price dynamics are given by the following: 

dS't = ^/yt + a^in StdW^ and dy^ = {a- by^)dt -b dZ^, y^ > 0. (40) 

For any a and 6 > 0, the last equation admits a unique strong solution y^ > 0. Note that the 
Lipschitz condition in Theorem 5.1 is violated, but by using the stopping times technique, we can 
directly show that Ey* < oo. Hence, this implies that condition (Ci) is satisfied for model (40). 

Similarly, we can check that (Ci) also holds for Ball-Roma models [3] or, more generally, for a 
class of processes with bounded diffusion satisfying the following condition. 

(A) There exist positive constants a, b, and M such that 

yFi(t,y) < a — by"^ and \F 2 (t,y)\<M, for all t>0,ySK. 

Proposition 5.1. Under condition (A), there exists a > 0 such that Ee“*^ < oo. 

Proof. The proof uses the same method as in Proposition 1.1.2 in [20]. □ 

Scott models: Suppose that volatility follows an Orstein-Uhlenbeck, as in Stein-Stein models, 
and the function a takes the exponential form 

d5't = (e'^*'‘ -b a^iJStdWy'> and dy^ = (a - foyjdt + dZ^, (41) 

where a,b and > 0 are constants. Here, 5 > 0 is chosen such that 25 < o, defined as in 
Proposition 5.1. Clearly, a{y) = e^^ + crand condition (Ci) is fulfilled because 

E sup |CT(y)p<2a^.„ + 2E(e2^1{| l>i}) <oo. 

0<t<l 

Numerical result for the Hull-White model: We provide a numerical example for the approx¬ 
imation result of Lepinette’s strategy in the Hull-White model (37). By Theorem 3.3, the corrected 
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n 

gain/loss 

corrected error 

lower bound 

upper bound 

price 

strategy 

10 

0.1523845 

-0.2225988 

-0.2363122 

-0.2088854 

0.7914033 

0.9013901 

50 

0.2966983 

-0.0596194 

-0.0670452 

-0.0521936 

0.9399330 

0.9706068 

100 

0.3086120 

-0.0288526 

-0.0350141 

-0.0226911 

0.9746527 

0.9875094 

500 

0.2955755 

0.0032387 

-0.0005821 

0.0070594 

0.9991733 

0.9995891 

1000 

0.2851002 

0.0012409 

-0.0021596 

0.0046415 

0.9999300 

0.9999652 


Table 1: Convergence for Lepinette’s strategy with = 0.01, g = 2. 


n 

gain/loss 

corrected error 

lower bound 

upper bound 

price 

strategy 

10 

0.2859197 

-0.0744180 

-0.0813544 

-0.0674816 

0.9246420 

0.9659700 

50 

0.3172523 

-0.0069238 

-0.0115426 

-0.0023049 

0.9921661 

0.9962377 

100 

0.3033519 

0.0007474 

-0.0030916 

0.0045864 

0.9984346 

0.9992385 

500 

0.3618707 

0.0001296 

-0.0024741 

0.0027333 

0.9999977 

0.9999989 

1000 

0.3334375 

0.0003996 

-0.0020559 

0.0028550 

1 

1 


Table 2: Convergence for Lepinette’s strategy with = 0.001, p = 4. 


replication error is given by F" — niax(S'i — K,0) — 77 min(S'i, K), where ry = 1 — K^(7{yi)g~^^/8/n. 

The difference F/* — max(S'i — K, 0) can be seen as the gain/loss of strategy y". For a numerical 
evaluation, we simulate N = 500 trajectories in a crude Monte-Carlo method, where the correlation 
coefficient of the two Brownian motions is 0.05 and the other initial values are Sq = K = 1, = 2, 

CTinin = 2,0 = —2 and 6=1. For each value of n, we estimate the average value of the corrected 
error and give the corresponding 95% intervals defined by lower and upper bounds. Initial num¬ 
bers of shares held are given in the last column of Tables 1 and 2. It turns out that strategy 7 " 
converges to the buy-and-hold strategy and the option prices approach the super-hedging price 
Sq. We also see that the convergence of the corrected replication error to zero is somehow slow. 
In fact, increasing values of g can provide a faster convergence, but this unexpectedly leads to 
super-replication more rapidly. 

We now provide a numerical illustration for the quantile hedging result of Proposition 4.2. For 
simplicity, suppose that a{y) = sin^( 2 /) -I- 0.1 and that y follows a geometric Brownian motion as 
above. To compare the reduction factor 1 — 5^ with powers of significance level e, we compute 
(1 — 5e)e~'' for 0.001 < £ < 0.1 and 0 < r < 0.1, with k * = 0.001. Then, (31) is confirmed by 
the simulation result (see Figure 2a). The simulation also shows that the option price inclusive of 
transaction costs is 1 — 0.385 = 0.615, which is cheaper than the super-hedging price = 1, for a 
shortfall probability less than 0,1%. Of course, it is reasonable to replace Sq by the option price 
inclusive of transaction costs C(0, Sq). The simulated reduction in the option price (1 —(5e)C(0, Sq) 
is then given in Figure 2b. 


6 High-frequency markets 

We now assume that purchases of the risky asset are carried out at a higher ask price S^ + e^, 
whereas sales earn a lower bid price S^ — e^. Here the mid-price S^ is given as in model (19) and 
Ej is the halfwidth of the bid-ask spread. Then, for any trading strategy of finite variation ip^, the 
wealth process can be determined by 


Ft = Fo + 






(42) 
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(a) Reduction factor 1 — 5^ and powers of e 


(b) Reduction amount and e 


where |'0| is the total variation of il’t- Observe that the first two terms are the classic components 
in frictionless frameworks, and respectively describe the initial capital and gains from trading. The 
last integral in (42) accounts for transaction costs incurred from the trading activities by weighting 
the total variation^ of the strategy with the halfwidth of the spread. 

For optimal investment and consumption with small transaction costs [21], the additional terms 
should be added in the formulation of V). In such cases, approximate solutions are usually deter¬ 
mined through an asymptotic expansion around zero of the halfwidth spread £, where the leading 
corrections are obtained by collecting the inputs from the frictionless problem. 

In this section, we are only interested in replication using discrete strategies in Leland’s spirit. 
Assume that for replication, the option seller applies a discrete hedging strategy , revised at 
n dates defined by ti = g{i/n) as in Section 3. The corresponding wealth process is now given by 


Vn,e ^ yn,e ^ 


f ^ et,\-tptf - I ■ 

“'O i=i 


(43) 


To treat the risk of transaction costs, we again apply the increasing volatility principle. Note that 
in high frequency markets, the bid-ask spread is, in general, of the same order of magnitude as 
price jumps^ . Hence, should be of the form for some positive constant Then, 

this case corresponds to the Leland-Lott framework with a = 1/2. 

In our context, the method in Section 3 is still helpful when replaced by Leland’s or 

Lepinette’s strategy. 

Proposition 6.1. Let and assume that the adjusted volatility is of the form 

o-s in (20). For both Leland’s and Lepinette’s strategies, the sequence of replicating 
portfolio values converges in probability to h(S']^)-|-min(S']^, AT) = S'!. In particular, — 

Si) converges to a mixed Gaussian variable as n —>■ oo. 

Proof. The proof is a direct consequence of Theorem 3.1, because the total transaction cost now 
converges to zero. □ 

Note that the case a = 0 studied in Section 3 corresponds to the assumption = k^S^, for 
some constant k* . This specific form means that the market is more illiquid and the bid-ask spread 


^It is important to know that the classical Black-Scholes strategy is not finite variation. 

®We thank an anonymous referee for pointing out the correspondence of the case a = 1/2 to this setting. 
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is now proportional to the spot price in every trade. Therefore, approximate hedging results for 
this case are the same as those in Section 3. 

We conclude the section by supposing that the stock spreads remain constant at all times, 
regardless of the current stock price. In other words, for some positive constant k*. 

Intuitively, transaction costs are now based on the volume of traded shares, instead of the traded 
amount of money as in the literature and Section 3. It is interesting to see that our methodology 
still works in this case. The following result is just an analog of Theorem 3.1, with a small 
modification to the limiting value of transaction costs, defined by 


Jo{x,y,Q) 





a{y)g 


\n{x/K) 

2A 



(44) 


where Z ^ A/’(0,1) independent of Si,yi 

Proposition 6.2. Suppose that > 0 and = Q\/nf'{t). For Leland’s strategy under 

condition (Ci), the sequence {V^’^ — h{Si) — min(S'j^, if) + k* Jo(5'i, 2 /ie)) weakly converges to a 

centered mixed Gaussian variable asn ^ oo. Furthermore, for Lepinette’s strategy, {yi ~ ~ ~ ^o) min(S']^,. 

weakly converges to a centered mixed Gaussian variable, where rjo = cr(2/i)£'“^5'f ^ a/S/tt. 

Proof. The proof is similar to that of Theorem 3.1 (see Section 7). □ 

Remark 11. When p —> oo under condition (C 2 ), one obtains an improved-rate version of Propo¬ 
sition 6.2, as in Theorem 3.2. 


7 Proofs 

Our main results are proved in the following generic procedure. 

Step T. Determine the principal term of the hedging error. In particular, we will show that the 
gamma term converges to 2min(S'i, if), while the cummulative transaction cost approaches J 
defined in (24). Both convergences are at rate 

Step 2: Represent the residual terms as discrete martingales. To this end, stochastic integrals will 
be discretized by following a special procedure set up in Subsection 7.2. 

Step 3: Determine the limit distribution of the normalized replication error by applying Theorem 
7.1. This result is the key tool, but we need in fact special versions adapted to our context. These 
will be explicitly constructed in Subsection 7.3. 


7.1 Preliminary 

Note first that C{t,x) and its derivatives can be represented as functions of and x, where 

Aj = Ag(l — t)^ := AQi/(t) and \=fLQ\/ri. (45) 


Moreover, the function tp{X,x), which appears in all fc-th (fc > 2) degree derivatives of C with 
respect to x and derivatives in time via the relation ( 8 ), is exponentially decreasing to zero when A 
tends to zero or infinity. This property motives our analysis in terms of variable A. In particular, 
let us fix two functions h, I* and let 1 < mi < m 2 < n be two integers such that G = \v{g{m 2 /n)) 
and I* = \v{g{rni/n)). Then, all terms corresponding to index j ^ [mi,m 2 ] can be ignored at 
a certain order which depends on the choice of G and I*. For our purpose, the desired order is 
9^ ~ Aq^. Therefore, we take, for example, G = 1/ In^ n, I* = In^ n and define 


m, = n — 


n(r/Ao)^/^^+^^ 


+ 1 and 


rrio = n — 




(46) 
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where the notation [cc] stands for the integer part of a real number x. Below, we focus on the 
subsequence {tj) of trading times and the corresponding sequence (A^), defined as 

= 1 - (1 - j/nY and = Ao(l - , rn^<j< m^. (47) 

Note that (f^) is an increasing sequence taking values in where = 1 — = 

1 — {I*/XqY^, whereas (A^) is decreasing in [Z*, T], Therefore, we use the notations Atj = tj —tj_i 
and AAj = \j-i — Aj, for mi < j < m 2 , to avoid recopying the negative sign in discrete sums. 

Below, Ito stochastic integrals will be discretized through the following sequences of independent 
normal random variables 


- W, 


^10 = 


( 1 ) 


\/^i 


and Z 2 j = 


^(2) _ p^(2) 

j _ j —1 

“ ^i-1 


(48) 


We set 


p{\x,y) 


Q fhi{x/K) 
a{y) V 2A 


(49) 


and write for short Pj-i = p{Xj_i,S^ ^,yt. j^)- This style of reduced notation is abusively applied 
for functions appearing in the approximation procedure. Define 


^3,j — l^i.i +Pj-il “E {\Z^ j +Pj-i| \T j_i), 

^4,^. = \z^Y - E {\z,Y I = \z,Y - v^. 


(50) 


The sequences (Z^ j) and (Z^^) will help to find the Dood decomposition of our approximation 
terms. In order to represent the limit of transaction costs, we introduce 


G{a) = E (|Z + a|) = 2p{a) + a (2$(a) - 1), 

A(a) = E (|Z + a| - E |Z + a\f = 1 + Y- G^{a), 


(51) 


for a S K and Z ^ JV{0, 1). We also write o{aY^) for generic sequences of random variables (Ar„) 
satisfying P - lim„_,^ = 0. 
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7.2 Approximation for stochastic integrals 

For any L > 0, we consider the stopping time 

T* =tI= inf {t > 0 : a{y^) + W{yt)\ > L} , (52) 

and denote by S* = and y^ = y^^^t the corresponding stopped processes. We provide an 

approximation procedure for ltd stochastic integrals through the sequences {Z^ j) and {Z 2 j). The 
discrete approximation concerns the class of functions satisfying the technical condition, 

(H) A : IR_|_ X x M —> K is a continuously differentiable function satisfying the following: there 
exist 7 > 0 and a positive function U such that for any x > 0, y G K, 

supmin(A'’', 1)|A(A, X, y)| < 17(x, y) and sup E y*) < 00, 

A>0 o<t<i 

for any —00 < m < +00, r > 0 and L > 0. 

Remark 12. We can check directly that for k > 2, d^C{X,x) = 

where P is some polynomial. Therefore, all functions Aq appearing in the approximation below are 
of the form X~^/‘^x"^a{y)P(hi{x/K)), where a can be a power of a or of its two first derivatives 
a' , a" . In bounded volatility settings, it can be shown with some computational effort (see e.g., 
[9, 24 , 27]) that 

sup ES”™ In^'^ St < 00, for any m G K, r > 0. (53) 

0 <t<l 

The latter property is, however, not always fulfilled for SV models with unbounded volatility. In 
fact, it has been demonstrated in [2, 28] that the stock price does not admit integrable moments in 
general SV markets, unless some natural conditions are imposed on the correlation and the volatility 
dynamics coefficients. Thus, asymptotic analysis using if estimates as in the existing works may be 
impossible in general SV frameworks. Nevertheless, note that (53) is true for processes stopped by 
T*. Below, the approximation analysis will be established in the sense of convergence in probability, 
in order to avoid this integrability issue. 

For simplicity, we use the notation S = {S,y). The following technique is frequently applied in 
our asymptotic analysis. 

Proposition 7.1. Let A{X,x,y) = Ao{X,x,y)(p{X,x), where Aq = Ao{X,x,y) is a function satis¬ 
fying (H). Then, for i = 1,2, 

1 \ m2 

A(A„ ^JdW« j dt = ^ A^_, Z.^^AX^ + o{e-^), (54) 

where 9^ = , A^ = A{X^,S^ ) and A{X,x,y) = A{z,x,y)dz. 

Proof. By making use of the stochastic Fubini theorem, we get 

4 = £ 4" Al(A„ dt = Al(A„ 5Jdt) dW« . 

Then, changing the variables v = X^ for the inner integral yields 

A{\, SJdt = A{v, 5Jdu = :4(A„, SJ - A(Ao, SJ. 
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In other words, /„ = /i „ - where ^ i„ = A(A„,5„) and = 

^(Aq, 5^) dlV^®). Moreover, we have 

A,n = f + r := + i?2,„ + i?3,„ • (55) 

•to -tt* -tt. 

Let e > 0 and b > 0. One observes that P( 6 *„|/ 2 ,nl > ■s) is bounded by P{t^ < 1) + P(^nl-^ 2 .nl > 
e, = 1). By condition (C^), we have 

limsupP(T* < 1) = 0. (56) 

L—¥00 

In view of (H), one deduces Al(Ag,x,j/)| < Cy/KU{x,y)e~^o/^, where U{x,y) = x~^/‘^U{x,y). 
Now, putting Al = and /*„ = A* dWb)^ one has P( 6 »„|/ 2 ,„| > e, r* = 1) = P( 6 »„|/*„| > 

e). Using the Chebychev inequality, we obtain 

P( 0 J/ 2 „I > £) < sup EU^iS;). 

’ ’ 0<t<l 


Hence, due to condition (H), l 2 n — as n —>■ oo. Similarly, taking into account that 

< Ajj < Aq for 0 < u < U, we get „ = o(6~^). 

Next, let us show the same behavior for the last term in (55). Indeed, for some fixed rj > 0 and 
L > 0, one has 

P (0Ji?3,nl > e) < P {dn\R3,n\ > S, + P , (57) 

where L^ ^ ^ = {inft^<u<i | ln(S'„/7f)| > r;, = l}. Then, by taking into account Lemma A.3 

and the integrability condition (Ci), one gets 

lim^^olim„^^limL^ooP(r)= ,^_^) = 0. 

On have A = A* and 

pOO 

ia:i < u{s*j / ( 1 + s*jdz < u{s:)f:, 


where /„* = + z-^)e-^^R^^'>-^/^dz. Set Lg,^- = {|ij < UiS^Jf:}, = A^l, 

and i?3 A*dW^®i. By Chebychev’s inequality, we obtain 


P( 0 „|i? 3 ,J <02e-2 / EiA*Jdu<ele-^ sup EU^iS*J {f:)^du, 

Jt 0<ii<l Jt 


which converges to zero as (/*)^du < Hence, i?3_„ = o{ 9 ~^). It remains to discretize 

the integral term i?2 « the sequence (Zij). The key steps for this aim are the followings. First, 
we represent i?2,„ = A„dIU^*) = ^ A^dlU^*^ and replace the Ito integral in the last 

sum with Aj_iZijy/Atj. Next, Lemma A.l enables us to substitute y/Atj = p“^AAj into the last 
sum to obtain the martingale Mm2 defined by Mk = Q~^ X)j=m Aj-iZijAXj. We need to show 
that |i?2,„ - MmJ = o(6»-i) or equivalently, Ejlm, Bj,n = o{ 9 Y): where A^^^dlUd) 

and A^ ■ = A(A„, S^) — A(A ,3, .Sj ). For this aim, we introduce the set 


F 


2,b — 


sup sup(|A(z,5„)| + \d^A{z,SJ\ + \dyA{z,SJ\) 

t* <U< 1 
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Then, for any e > 0, P J2T=m.^ ^j,n\ > is bounded by P(r 2 j,) + P(r* < 1) + zun, where 
ru„ = P (ej -Bj.nl > £, ^2,b, T* = iy Let where 

= ^«jl{|T„,d<KIA„-A,_d+|s;-s; l+lvl-y; DL 

Then, ti7„ = P ^i,nl > e). '^'hich is smaller than £"^6*^ E by Chebychev’s 

inequality. Clearly, E i?| ^ is bounded by 


36^ 


((A„ - A,_,)2 + E(5: - s: , + E« - yl ^ )^)du < {AX,f + {At,f 


up to a multiple constant. Consequently, 0'^Y^=m^ (^A^)^ + (At^)^, which 

converges to 0 by Lemma A.l and condition (C 2 ). On the other hand, by Lemma A.4, we get 
lim^_,,,,^ lim„_,,gQP(r 2 j) = 0. The proof is complete. □ 


7.3 Limit theorem for approximations 

We first recall the following result in [15], which is useful for studying asymptotic distributions of 
discrete martingales. 

Theorem 7.1. [Theorem 3.2 and Corollary 3.1, p.58 in [15]] Let A4„ = ® zero-mean, 

square integrable martingale and c be an a.s. finite random variable. Assume that the following 
convergences are satisfied in probability: 

n n 

-^0 for any S >0 and ^ ^. 

Then, (Al„) converges in law to X whose characteristic function is Eexp(— i.e., X has a 
Gaussian mixture distribution. 

In this subsection, we establish special versions of Theorem 7.1. In fact, our aim is to study the 
asymptotic distribution of discrete martingales resulting from approximation (54) in Proposition 
7.1. First, we define 

k 

A4fe = ^ Vp < fc < m 2 , (58) 

j=m^ 

where A.^j_i Z^jAXp ^ = A^{Xj,Sf, ) and Zij defined as in (48) and in (50). To 

describe the asymptotic variance of (M), we introduce the following function 

L(A,a;,y) = Al{X,x,y)-\-2A^{X,x,y)A^{X,x,y){2<^{p) - 1) 

+ Al{X,x,y)A{p)Al{X,x,y), (59) 

where p is defined in (49). Set 

+ and /i= (/r-l)/(/i + l). ( 60 ) 

Proposition 7.2. Let A° = A°(A, x,y), i = 1,2,3 be functions having property (H) and A^{X, x, y) = 
A^{X,x,y)ip{X,x). Then, for any fixed g > 0, the sequenee {n^M.rn 2 )n>i weakly converges to a 

mixed Gaussian variable with mean zero and variance defined as = <r^(<Si) = figt^ A^L(A, .S]^)dA. 
The same property still holds if some (or all) of the functions Ai areoftheform A[(z, x,y)(p(z, x)dz. 
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Proof. Note that the square integrability property is not guaranteed for (vj). To overcome this 
issue, we consider the “stopped version” {v*), which is obtained by substituting S^._^ by Sf ^ 
in A^, i.e., v* = S*_) Let Ml = the corresponding stopped 

martingale. First, we show, throughout Theorem 7.1, that for any L > 0, this martingale weakly 
converges to a mixed Gaussian variable with mean zero and variance To this end, 

setting Ti = {inft.<„<i | \ii{Sl/K)\ > rj} and a* = E |we obtain 


^ 2/3 I 


j=m^ 


> £ < P 




E r,, 1 +P(rj^^). 




(61) 


It suffices to show that the first probability on the right side of (61) converges to zero. Indeed, 
from the proof of Proposition 7.1, one observes that on the set 


max 
2 = 1 , 2 ,3 


U.(A„,5:)|<C/(5*)(1 + Ay), 




(62) 


for some 7 > 0 and U{S) = S ^/^U{S). Set v* = z;*lp and a* = E (u*^l/|o*|> 5 \|.F _^), where 

Tg^^. = I max^ |A,(A„,^:)| < C/(5:)(l + A;^)|. 


1 < 2<3 


We have 


P a*| > e, T^^^ ^ E ^ ^ E 


^=mi 


j=mi 




by Markov’s inequality. By using Chebychev’s inequality and then again Markov’s inequality, we 
observe that 

Ea; < y^yP(|u;| ><5) < S-^Ed*/ < 96-^{l + X-^fiAX^fEU^Sl) Y 


Note that has bounded moments. Then, by using (62), we obtain 

bounded by A((^)"^(^Aj)'^, which converges to zero by Lemma A.l. 

We now verify the limit of the sum of conditional variances E(u*^| _;^). Set v* . = A* Z^ AA -. 

3 3 '^^3 ^t 3 i- '3 J 

As Z'jj and ^27 are independent, E (^l -v^ .\F ^= E -v^ = 0. It follows that 


= E(q^,-|^,_i) + E(u2*^^.|J-,_i) + E(q;^^.|J-^._i) + 2E{vlY2j3^3-i)- 


Now, observe that for Z ~ A(0,1) and some constant a, E(Z \Z + a|) = 2$(a) —1 and E (Z + a)^ — 
(E|Z + a|)^ = A(a). On the other hand, AA^ = n“^^(l + o(l))/i by Lemma A.l. There¬ 

fore, 

n^^E{vf\F^_,) = (1 + o(l))Ae^ XY 

By Lemma A. 5, Y^3^ E{v*^\F^_i) converges in probability to c*^(L). Thus, weakly 

converges to ^(0, <^*^(L)) by Theorem 7.1. Moreover, property (56) implies that for any i5 > 0, 

lim lim p(n^|M^^- m;;; I > 5) = 0 . 

L—¥oo n—^oo 
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Therefore, by taking into account that converges a.s. to as L —)■ oo, we conclude that 

converges in law to A/’(0,^^). This completes the proof. □ 

Next, we study the asymptotic property of the following martingale 

k 

M-k = ^Ij + ^2,j-l ■^2,j + ^4,i) (63) 

i="ii 

The limiting variance will be defined throughout the function 

L(A,a;,y) = A\{\x,y) + Al{\x,y) + (1 - 2/Tv)Al{X,x,y). (64) 

The following result is similar to Proposition 7.2. 

Proposition 7.3. Let = 7l°(A, x,y), i = 1, 2,4 be functions having property (H) and A^{X, x, y) = 
A^{X,x,y)ip{X,x). Then, for any fixed g > 0 the sequence M.m^)n>i weakly converges to a 

mixed Gaussian variable with mean zero and variance given byg^=jl L(A, 5]^)dA. 

The same property still holds if some (or all) Ai are of the form A^(z,x,y)(p(z,x)dz. 

Proof. The conclusion follows directly from the proof of Proposition 7.2 and the observation that 
E = E(|Zj^ ^1 — A/ 2 / 7 r)^ = 1 — 2 / 7 r, and E ^^4 ^■) = 0, for i = 1, 2 and < j < ^ 3 . □ 

In the rest of the subsection, we establish a limit theorem for a martingale of the following form 

k k 

^k= + ^3j-l ^3,j) := < fc < TO2, 

where A^{X,x,y) = A){X, x,y)ip{X, x) and A), i = 1,3 are functions having property (H). The 
following result is helpful for the case when g diverges to infinity as in Theorem 3.2. 

Proposition 7.4. Under condition (C 2 ), the sequence weakly converges to a 

mixed Gaussian variable with mean zero and variance = fi X^ L(A, 5 ' 4 )dA, where L(A, x, y) = 
A‘){X, X, y) + 2 A 14 (A, X, y)A^{X, x, y) + ^^(A, x, y). The same property still holds if some (or all) Ai 
are of the form A)(z, x,y)^(z, x)dz. 

Proof. We determine the limit of conditional variances of gt^ observe that 

= A(1 + 0(1)) Af_4 Q(A^._i, JAA^-, (65) 

where Cl{X,x,y) = Al{X,x,y) + Al[X,x,y) K{p) + 2A4(A, a;, y) 2 l 3 (A, a;, j/) (24>(|p|) - 1). Moreover, 
it can be checked directly that the function G(-) defined in (51) satisfies the following inequalities: 
|a| < G{a) < |a| + 2ip{a) , for any a S K. This implies that |A(a) — 1| < 4|a|95(a) + 
hence, sup^^j^ |A(a)| < 00 . Note also that Q —>• L a.s. as n —>■ 00 because p{X,x,y) —>■ 00 as 
g = g{n) —>■ 00 , for any x > 0 and A ^ 21n(x/Ar). Using Lemma A.5, we claim that the sum of 
the terms on the right-hand side of (65) converges in probability to The proof is completed by 
running again the argument in the proof of Proposition 7.2. □ 

7.4 Proof of Theorem 3.1 

We first observe that 14 „ approaches 2 min( 5 ' 4 ,Ar) at order In particular, setting 14 „ = 
(Sf(p(Xf, S'J — S' 4 ^(A(, 54 )) dt and changing variables v = t?sds, we can represent 
li,n = *^1 Iq° Si) dx -I- 14 „ -b o{0~^). The hrst integral in the right side converges a.s. to 
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2min(5'i,iir) by ( 12 ), while„ is approximated by dt, where 

iJ = (2“^A“^/^ — ln(a;/it'))^(A, cc). The discretization technique of Proposition 7.1 can be 

applied to replace the latter double integral by defined by 

k 

^i,fc = mi <k <1712, (66) 

j^mi 

where H{X,x) = 12 — ln{x/K))(p{z, x)dz. We summarize the asymptotic form of 

in the following. 

Proposition 7.5. If g either is constant or satisfies condition (C 2 ) then, 

P- lim 9^\f -2mhi{S^,K)-U^ \=d. 

n—yoo 


Next, we claim that /2 ^ = o{9^^). 

Proposition 7.6. If g either is a positive constant or satisfies condition (C 2 ), then 0„72 n 
verges to zero in probability as n —>■ 00 . 

Proof. See the Appendix B. □ 

Let us study the trading volume J„. First, it is easy to check that for any ?; > 0, 1 — <i>(n) < 


^g:>{v). Now, observe that I 7 " — 7 " 


< |1 - 7 "| + |1 - 


which almost surely con¬ 


verges to zero more rapidly than any power of n when infj^<j<;„Aj > I* i < mi. The 

same property can be deduced for the case sup. A < /* <;=^ i > m 2 . To see this, we note that 
for Xu < I*, Su{uj) = converges to 5'i(a;) as n —)■ 00 uniformly in A„ G [0,/*], 

for any tv outside the zero probability set {S'! = K}. Therefore, one can truncate the sum 
and keep only the part corresponding to index mi < j < m 2 . In other words, Jn is ap¬ 
proximated by Ji „ = ZTiru, ^t. |A$,| . Putting &, = |A$,| — (Pj_i l^v l, we can represent 


A.n = Jin + ^i,n + £ 2 ,™ : where ' 

£ 2 ,™ = JiTirm by In view of (70) and condition (C 2 ), ei^„ = o{0-^) as 


AS, 


_ y^m2 

“ Aij=m, -j_i 

n —>■ 00 . Furthermore 


Aj$| and 


by using the Taylor expansion, we obtain | £2 „ | < *S'sup Yiy-rn-, 
where = supp^^^j^ S,. Now, by taking into account that 




A-,! 


up to a multiple constant. 


^ + (AA-Af)N(AJ_7-A7'/^)‘ 


up to a multiple constant and using condition (C 2 ) together with (70), we get 1^2 nl = 
Next, by using Ito’s Lemma and the substitution \j = Aq(1 — we replace ^ by 


Jo 'n - 


m 2 m 2 


■ 'J 

j^mi 


^j-i • 


j=mi 


Cy K^= g ^a{y,^_^)Zi^^+q^_i, 


(67) 


where q is defined in (10). We need to determine the limit of Jn throughout the Doob’s decompo¬ 
sition of Jo „ w.r.t. the hltration (Pf) . To this end, note that 

A J1 mi<j<m2 

E(<;,I7_,) = A-y"7_y,_,AA, E(|x,||J)_,), 

where E(|>r^| |J(,_i) = g~^a{y, ^)G{pj_i) Dj_i and G{p) defined in (51). Let 

m 2 

B{X,x,y) = X~^^^xip{X,x)D{X,x,y) and B^_iAXj. ( 68 ) 


j^mi 
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We observe that J 2 „ = J 3 „ + 1^2,m 2 ^ 


where 


^ 2 .fe = and := \xj\- D^_^. 


(69) 


By substituting by Si everywhere in we write J 3 „ = J 4 n ^5 m where J^n — 

i3(A,_i,^i)AA,-, h,n = S;-iAA- and B*_^ = B(A,_i,5,’_J-i?(A,_3,5i). Then, 

by Lemma A.2, we can check that J 4 „ converges a.s. to g) at rate 6 *„. Now, an application 

of Ito’s Lemma for yields stochastic integrals with respect to the Wiener processes. Owing 
to Proposition 7.1, the sum of these integrals can be approximated by defined by 

2 k 

^3,k = XI X mi<k< m 2 , 


where Qi = (xa{y)dxB + rF 2 {t{X),y)dyB)dz and Q 2 = Vl — r^F 2 (t{X),y) /“ dyBdz and 
t{X) = 1 — (A/Ag)"^^. The asymptotic form of J„ is summarized in the following. 

Proposition 7.7. For any fixed g > 0, 

P- lim I J„ - 7(53,2/3,0) - (Z^2,m2 +^ 3 .mJI = 0 - 

n—yoo 

Now, the martingale part of the hedging error is given by defined by 


1 ^ ^ ^ 

^k = 2^1,k - «^*(^2.fc +^ 3 .fc) = X mi<k< m2, 

j=mi i=l 

where A 3 = —(j{y)xH/2, A 2 = k^Q 2 and A 3 = —K^(j{y)X~^/‘^Xip{X,x). It is easy to see that 

the assumption of Proposition 7.2 is fulfilled for A^, i = 1,2,3. Hence, the sequence 

converges in law to a mixed Gaussian variable by Proposition 7.2, which proves Theorem 3.1. □ 


7.5 Proof of Theorem 3.2 

When g —>■ 00 under condition (C 2 ), the approximation for J„ is slightly different. In particular, 
observe first that for any 6 G K, E jaZ + b\ can be approximated by 6(2<I>(6/a) — 1) as a —>■ 0. 
Therefore, we can replace J 3 „ in ( 68 ) by the sum ^ where B{X,x) = 

X~^^^ xip{X,x)q{X,x)^{gq{X,x)), with ^{q) = 2^{gq) — 1 and q{X,x) defined in (10). Puting 
An = %-i = - B{Xy_i,Si), we represent := 

n = EE B*_, AA . Now, using Lemma A. 2, we can directly show that | J 4 „ —J*( 53 )| = 
o{9~^). Furthermore, owing to Ito’s formula, we replace A-i A /t SxB{X^_i, Sfi)dS^. Direct 
calculations yield that 

d 2 ,B = A"^/^ fi{X,x)[-2q^{X,x)^{X,x) + ^$(A,x) + |( 0 ( 0 g(A, x))]. 

Clearly, ^{gq) —)■ sign.{q) and (p{gq) —)■ 0 as p —/■ 00 . Now, using Proposition 7.1, we can approx¬ 
imate J 5 „ by defined by ^ E^Xi N{X,x) = 

x) {-2q^{z, x) + l/ifiz)) sign(g(z, x))dz. The asymptotic representation of the trad¬ 
ing volume is summarized in the following. 
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Proposition 7.8. Under conditions (Ci) — (C 2 ), 

P- lim 0J4- 

n—yoo 

Now, the martingale part of the hedging error is determined by 

k 

■^k = ~ +^3,k) = ^ + ^3,j-lZ3j)AXj, 

j=mi 

where Ai, i = 1,2 are explicitly determined and satisfy the assumption of Proposition 7.4. Then, 
in view of Theorem 3.2 is proved throughout Proposition 7.4. □ 

7.6 Proof of Theorem 3.3 

The key technique in Proposition 7.1 can be used to obtain a smart martingale approximation for 
the sum J2i>i ^5'*, S'„)du. 

Proposition 7.9. If g either is a positive constant or satisfies condition {C 2 ), t/ien I /2 m 2 1 = 
o{9~^), where Y{X,x) = K)ip{z,x)Az and 

k 

^i,k = 8~^ mi < fc < TO 2 . 

j^mi 

Proof. The proof follows from the substitution AS^. by g~^a{y^._^)S^._^AX^. as in Proposition 
7.1. □ 

Let us now study the trading volume J„ by following the procedure in the approximation of 
Jn- First, by Ito’s lemma, 

^ S'„)cr^(j/„)S'^dw, 

where the time-correction, which involves the term qj-i in the formula of Xj defined by (67), has 
been removed. We now approximate by 

m2 

Ji,n = 8~^ Bj_j^\Zj^j\AXj and B{X,x,y) = (j{y) x X~^^'^ip{X,x). 

j^mi 

As E|Z| = \/2l'K, for Z ^ A/’(0,1), the Dood’ decomposition of Ji „ is given by </2 n + ^2 m j 
where = g-^ ^_^AX^ and ^a.m^ = E”Lmi putting 

- ■B(Aj_i,.5i), we write = ^A.n^^ 3 ,-a^ where 

m 2 m 2 

74,„ = ^3,n = 51 

3—rax i=mi 

Observe that J 4 ^ converges a.s. to 77 min(S'i, iF) by Lemma A.2 and (12). We now find the suitable 
martingale approximation for J 3 By Ito’s formula, B^i can be replaced by It Qi(Aj_i, <5„)dW'i*\ 
where = a{y)xdxB + rF 2 {t{X),y)dyB and = Vl — r^F 2 (t(X),y)dyB. Direct calculations 
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show that dxB = cr(y)(2 — A \n{Xl K))ip{X,x) and dyB = a'{y)\ Now, 

Proposition 7.1 can be applied to approximate J3 „ by the martingale 773^^2 > defined by 

k 

7^3,fc = Q ^ + ^2,j-i^2,i)AAj, mi <k< m 2 , 

j=mi 

for explicit functions A^,i = 1,2. The final asymptotic form of is given below. 

Proposition 7.10. If g is a positive constant independent of n then, 

P- lim 0 „|J„-? 7min(5 'i,X)-( 7/2 + 7^3 )|=0. 

n—¥oo 

Hence, the martingale part of the hedging error for Lepinette’s strategy is determined by = 

^i.mj +^1.7712 ~ + The latter martingale sum can be represented in the form 

k 

Adj. = Q ^ {Aij_iZi^j -\- Aji^jZjij_i + A2j_iZ2j)/X.\j, mi <k< m 2 , 

j=mi 

for explicit functions holding the assumption of Proposition 7.3. Then, >3 converges 

in law to a mixed Gaussian variable, which completes the proof. □ 


8 Conclusion 

We studied the problem of approximate option replication in SV settings using a new specification 
for adjusted volatility. Although our model employed a simpler adjusted volatility than in the 
previous literature, we obtain the same asymptotic results for both Leland’ and Lepinette’s strate¬ 
gies in general SV markets. A possible connection to high frequency markets with proportional 
transaction costs was also discussed. As an application, we showed that the option price inclusive 
of transaction costs can be reduced by adapting the theory of quantile hedging. Note that our 
approach is still helpful for more general settings, for example, when the friction rule admits a 
separate-variable representation [31]. This generalization includes the case where trading costs 
are based on the physical number of traded shares. Lastly, in a companion paper, we extended 
the method to multidimensional frameworks for European options with general payoffs written on 
several assets [32]. 
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Appendix 

A Auxiliary Lemmas 


Lemma A.l. There exist two positive constants Ci,C 2 such that 

Cin~^^vii{U) < inf jAAjj< sup \AXj\ < C 2 n~^^gi^ 

^1—i—^2 

where voix) = ^ Moreover, for any mi < j < m 2 , 

AXj = n~'^^g^^ iyo{Xj-i){l + o{l)) and AA^-= g(l-|-o(l)). 


(70) 


(71) 
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Proof. It follows directly from the relation (47). □ 

A technical condition (Hq): A : K_|_ —>■ IR js a continuously differentiable function having 
absolutely integrable derivative A' and 


lim 


|A(A)|dA + 


^+oo 


|A(A)|dA =0, 


where . 


The following result is straightforward to check. 


Lemma A.2. Let g either be a positive constant or satisfy condition (C 2 ). Then, for any function 
A satisfying condition (Hq), 


lim 

n—¥oo 


^2 .00 

Yl - / ^(A)dA 

J a 




= 0. 


(72) 


In particular, hm„_^^ 6^ 


ETJ^A(A^_M - J-A(A)dA 


= 0 . 


- J Jo 

Lemma A.3. For any K > 0, lim^_^g limsup^ ^ P(inf^<.^<^ | ln(5'„/A)| < e) = 0. 


Proof. It follows from the fact that conditioning on the cr-held generated by the volatility process, 
the log-price process In has Gaussian distribution. □ 

Lemma A.4. Suppose that Aq and its derivatives dxAf^, dyA^ verify condition (H). Set A(A, x, y) = 
Ao{A,x,y)lp{A,x), A{A,x,y) = /“ A{z ,x,y)dz and define 

r„ = sup {\d^A{z,r,d)\ + \d„,A{z,r,d)\ + \dyA{z,r,d)\) , 

{z,r,d)^ [Z^ ,l*]xl3 


where 3 ‘^max] ^ ymax] ‘^min '^max Vniin 

Vu, ymax = suP(.<„<t^ y„. Then, lim^^^ lim„^^P(r„ > 6) = 0. 


Proof. Let e > 0. On the set ^ = {inf(.<^<;^ | ln(S'„/A)| > e}, 

sup TiHA) < (27r)“^/^\/ Kr~^ exp{—e^/ (2q) — q/8}. 

Synin<r<S„,,,„ 

By condition (H), there exists 7 > 0 such that 

^ pOO ^ 

|A^(z,r,d)| < C|17(r,d)| / (g-^/^ + g 7 )e-eV( 29 )- 9 / 8 g^g < 

J Z 

where U is some function verifying supg^^^^ E,U{Sf) < 00 . For any 77 > 0 and N > 0, let 

r 2 ,, = { sup |G(r,d)-G(5i)| <r7}f]{|t7(5i)| < A}. 

{r,d)eB 


It is clear that |17(r, d)| < A -|- 7 on the set r 2 _^. Similarly, taking into account dxA{z,r,d) = 
—A{z, r, d), dyA{z, r, d) = dyAo{z, x, y)Ip{z, x) we deduce that both |i9;^A(z, r, d)\ and \dyA{z, r, d)\ 
are bounded on r 2 _,j by a constant independent of b. Now, for & > A -|- 77 -|- 20^^ ^, P(r„ > b) 
is bounded by 

P(rj J + P( sup |G(r, d) - U{Sl)\ > 77 ) + P{\U{Sl)\ > A) + P(r* < I). 

’ (r,d)eB 

By Lemma A.3, limg_,.g lim„_,, 3 Q P(r)^ g) —t 0. Thanks to the continuity of the functions St and j/^, 

one gets hm„_,,^ P (^sup^^^^^g |A(r,c?) — A(,S*)| > 77 ^ =0. Moreover, the integrability of U{S*) 

implies that PdA)^)")] > A) converges to zero as A —>■ 00 . By (56), P(r* < 1) converges to 0 as 
L —>■ 00 , which completes the proof. □ 
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Lemma A.5. Let A{X,x,y) = x,y)(p{z,x)dz, A = A^, where is a function having 

property (H). Then, for any 7 > 0, 


P — lim 

n—>-oo 


m 2 


j=mi 



A^A(A,5i)dA 


= 0 , 


where Sf = {Sf,yf). The same property still holds if A(X,x,y) = A^{X,x,y)(p{x,y) or is a product 
of these above kinds. 


Proof. We prove for the first case A{X,x,y) = A^{z,x,y)ip{z,x)dz, as the same argument 
can be made for the other cases. First, we split the expression under the absolute sign as 
AJ_iA(Aj_i,5i)AAj + Ejlmi where A^.„ = A(A^_i,5t._J - A(A^_i,5i) and 


A{X,x,y) = X^A{X,x,y). It is clear that for any {x,y), the function A{-,x,y) satisfies condition 
(Hq). Hence, ^(Aj-ij'5i)AAj converges a.s. to zero by Lemma A.2. It remains to show 

that P(|A„| > e) —)■ 0 for any given e > 0, but it can be done by the same way as in Lemma 


A.3. □ 


B Proof of Proposition 7.6 


The singularity of C at the maturity T = 1 requires a separate treatment. Let e„ = n g 
We then represent / 2 _„ = w^{t)dW^^'’zu„it)dW^^\ where zujt) = {jlf-C„,{t, St))a{y^)S^. 

Taking into account that I 7 " — C,^{t,S^)\ < I, we obtain lim„_^^ 0^ E vj^{t)dt = 0. Now 


put t = min(t-, I — e„). It then remains to prove that X]"-i ^(7^ ~ = o(0^^). 


>-1/2 


Let us introduce the discrete sums Wi{t) = -1 (At 

J — J- ^j — l^ J — X t 

^ 3(0 = where and x^ = \n{SjK). 

Clearly, I 7 ” — S'Jp < w.^{f) + W 2 {t) + w^{t). By taking into account that 


sup n sup 'E{x^ — xp f,j{t) < 00 and 

n,l<j<n 0<t<l 


sup < 00 , 

0<f<l 


we have 0^E rt;j^(t)dt < which converges to zero by (C 2 ). Now, the particular 

choice of e ensures that 0^E W 2 (t)dt < Cd^n~^X^^ , which tends to zero. The 
convergence for w^lt) can be shown in the same way. □ 


C Moments of Orstein-Uhlenbeck’s processes 

Lemma C.l. Suppose that cr(z) < 7(1 + |z|) for all z € K, for some constant 7 > 0. Let y^ he an 
Or stein-Uhlenheck process defined by dy^ = (a — byf)dt-t-dZf with some constants a and b > 0. Put 

= exp | 2 a 7 ^ l/sdsl and = b^(2j^(2b + a^)) ^. Then, EA^ < 00 for any 0 < a < 

Proof. Remark that {a — by)y < j{2b) — hy^ 12. Then, by adapting Proposition 1.1.5 in [20, p.24], 

we can show that Elj/jp™ < ml (2fb + a^fb^) for any integer m > 1. It follows that for any 
0 < a < a^, 

EA„ < ^ (a 27 ^)-(m!)-iE|y/™ < ^ (2/5 + (a27^)'" < 00 . 

771=0 771 = 0 

If y^ is mean-reverting then 5 takes very large values. Hence, it is possible to choose a > 3/2 -b ^/2 
as discussed in Remark 10. □ 
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